\(\int \frac {(A+B x) \sqrt {d+e x}}{(b x+c x^2)^2} \, dx\) [1242]
Optimal result
Integrand size = 26, antiderivative size = 158 \[
\int \frac {(A+B x) \sqrt {d+e x}}{\left (b x+c x^2\right )^2} \, dx=-\frac {(A b-(b B-2 A c) x) \sqrt {d+e x}}{b^2 \left (b x+c x^2\right )}-\frac {(2 b B d-4 A c d+A b e) \text {arctanh}\left (\frac {\sqrt {d+e x}}{\sqrt {d}}\right )}{b^3 \sqrt {d}}+\frac {\left (2 b B c d-4 A c^2 d-b^2 B e+3 A b c e\right ) \text {arctanh}\left (\frac {\sqrt {c} \sqrt {d+e x}}{\sqrt {c d-b e}}\right )}{b^3 \sqrt {c} \sqrt {c d-b e}}
\]
[Out]
-(A*b*e-4*A*c*d+2*B*b*d)*arctanh((e*x+d)^(1/2)/d^(1/2))/b^3/d^(1/2)+(3*A*b*c*e-4*A*c^2*d-B*b^2*e+2*B*b*c*d)*ar
ctanh(c^(1/2)*(e*x+d)^(1/2)/(-b*e+c*d)^(1/2))/b^3/c^(1/2)/(-b*e+c*d)^(1/2)-(A*b-(-2*A*c+B*b)*x)*(e*x+d)^(1/2)/
b^2/(c*x^2+b*x)
Rubi [A] (verified)
Time = 0.25 (sec) , antiderivative size = 158, normalized size of antiderivative = 1.00, number of
steps used = 5, number of rules used = 4, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.154, Rules used = {834, 840, 1180, 214}
\[
\int \frac {(A+B x) \sqrt {d+e x}}{\left (b x+c x^2\right )^2} \, dx=-\frac {\text {arctanh}\left (\frac {\sqrt {d+e x}}{\sqrt {d}}\right ) (A b e-4 A c d+2 b B d)}{b^3 \sqrt {d}}+\frac {\left (3 A b c e-4 A c^2 d+b^2 (-B) e+2 b B c d\right ) \text {arctanh}\left (\frac {\sqrt {c} \sqrt {d+e x}}{\sqrt {c d-b e}}\right )}{b^3 \sqrt {c} \sqrt {c d-b e}}-\frac {\sqrt {d+e x} (A b-x (b B-2 A c))}{b^2 \left (b x+c x^2\right )}
\]
[In]
Int[((A + B*x)*Sqrt[d + e*x])/(b*x + c*x^2)^2,x]
[Out]
-(((A*b - (b*B - 2*A*c)*x)*Sqrt[d + e*x])/(b^2*(b*x + c*x^2))) - ((2*b*B*d - 4*A*c*d + A*b*e)*ArcTanh[Sqrt[d +
e*x]/Sqrt[d]])/(b^3*Sqrt[d]) + ((2*b*B*c*d - 4*A*c^2*d - b^2*B*e + 3*A*b*c*e)*ArcTanh[(Sqrt[c]*Sqrt[d + e*x])
/Sqrt[c*d - b*e]])/(b^3*Sqrt[c]*Sqrt[c*d - b*e])
Rule 214
Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[-a/b, 2]/a)*ArcTanh[x/Rt[-a/b, 2]], x] /; FreeQ[{a, b},
x] && NegQ[a/b]
Rule 834
Int[((d_.) + (e_.)*(x_))^(m_)*((f_.) + (g_.)*(x_))*((a_.) + (b_.)*(x_) + (c_.)*(x_)^2)^(p_), x_Symbol] :> Simp
[(d + e*x)^m*(a + b*x + c*x^2)^(p + 1)*((f*b - 2*a*g + (2*c*f - b*g)*x)/((p + 1)*(b^2 - 4*a*c))), x] + Dist[1/
((p + 1)*(b^2 - 4*a*c)), Int[(d + e*x)^(m - 1)*(a + b*x + c*x^2)^(p + 1)*Simp[g*(2*a*e*m + b*d*(2*p + 3)) - f*
(b*e*m + 2*c*d*(2*p + 3)) - e*(2*c*f - b*g)*(m + 2*p + 3)*x, x], x], x] /; FreeQ[{a, b, c, d, e, f, g}, x] &&
NeQ[b^2 - 4*a*c, 0] && NeQ[c*d^2 - b*d*e + a*e^2, 0] && LtQ[p, -1] && GtQ[m, 0] && (IntegerQ[m] || IntegerQ[p]
|| IntegersQ[2*m, 2*p])
Rule 840
Int[((f_.) + (g_.)*(x_))/(Sqrt[(d_.) + (e_.)*(x_)]*((a_.) + (b_.)*(x_) + (c_.)*(x_)^2)), x_Symbol] :> Dist[2,
Subst[Int[(e*f - d*g + g*x^2)/(c*d^2 - b*d*e + a*e^2 - (2*c*d - b*e)*x^2 + c*x^4), x], x, Sqrt[d + e*x]], x] /
; FreeQ[{a, b, c, d, e, f, g}, x] && NeQ[b^2 - 4*a*c, 0] && NeQ[c*d^2 - b*d*e + a*e^2, 0]
Rule 1180
Int[((d_) + (e_.)*(x_)^2)/((a_) + (b_.)*(x_)^2 + (c_.)*(x_)^4), x_Symbol] :> With[{q = Rt[b^2 - 4*a*c, 2]}, Di
st[e/2 + (2*c*d - b*e)/(2*q), Int[1/(b/2 - q/2 + c*x^2), x], x] + Dist[e/2 - (2*c*d - b*e)/(2*q), Int[1/(b/2 +
q/2 + c*x^2), x], x]] /; FreeQ[{a, b, c, d, e}, x] && NeQ[b^2 - 4*a*c, 0] && NeQ[c*d^2 - a*e^2, 0] && PosQ[b^
2 - 4*a*c]
Rubi steps \begin{align*}
\text {integral}& = -\frac {(A b-(b B-2 A c) x) \sqrt {d+e x}}{b^2 \left (b x+c x^2\right )}-\frac {\int \frac {\frac {1}{2} (4 A c d-b (2 B d+A e))-\frac {1}{2} (b B-2 A c) e x}{\sqrt {d+e x} \left (b x+c x^2\right )} \, dx}{b^2} \\ & = -\frac {(A b-(b B-2 A c) x) \sqrt {d+e x}}{b^2 \left (b x+c x^2\right )}-\frac {2 \text {Subst}\left (\int \frac {\frac {1}{2} (b B-2 A c) d e+\frac {1}{2} e (4 A c d-b (2 B d+A e))-\frac {1}{2} (b B-2 A c) e x^2}{c d^2-b d e+(-2 c d+b e) x^2+c x^4} \, dx,x,\sqrt {d+e x}\right )}{b^2} \\ & = -\frac {(A b-(b B-2 A c) x) \sqrt {d+e x}}{b^2 \left (b x+c x^2\right )}+\frac {(c (2 b B d-4 A c d+A b e)) \text {Subst}\left (\int \frac {1}{-\frac {b e}{2}+\frac {1}{2} (-2 c d+b e)+c x^2} \, dx,x,\sqrt {d+e x}\right )}{b^3}+\frac {\left (2 \left (\frac {1}{4} (b B-2 A c) e+\frac {\frac {1}{2} (b B-2 A c) e (-2 c d+b e)+2 c \left (\frac {1}{2} (b B-2 A c) d e+\frac {1}{2} e (4 A c d-b (2 B d+A e))\right )}{2 b e}\right )\right ) \text {Subst}\left (\int \frac {1}{\frac {b e}{2}+\frac {1}{2} (-2 c d+b e)+c x^2} \, dx,x,\sqrt {d+e x}\right )}{b^2} \\ & = -\frac {(A b-(b B-2 A c) x) \sqrt {d+e x}}{b^2 \left (b x+c x^2\right )}-\frac {(2 b B d-4 A c d+A b e) \tanh ^{-1}\left (\frac {\sqrt {d+e x}}{\sqrt {d}}\right )}{b^3 \sqrt {d}}+\frac {\left (2 b B c d-4 A c^2 d-b^2 B e+3 A b c e\right ) \tanh ^{-1}\left (\frac {\sqrt {c} \sqrt {d+e x}}{\sqrt {c d-b e}}\right )}{b^3 \sqrt {c} \sqrt {c d-b e}} \\
\end{align*}
Mathematica [A] (verified)
Time = 0.55 (sec) , antiderivative size = 151, normalized size of antiderivative = 0.96
\[
\int \frac {(A+B x) \sqrt {d+e x}}{\left (b x+c x^2\right )^2} \, dx=\frac {\frac {b \sqrt {d+e x} (b B x-A (b+2 c x))}{x (b+c x)}+\frac {\left (4 A c^2 d+b^2 B e-b c (2 B d+3 A e)\right ) \arctan \left (\frac {\sqrt {c} \sqrt {d+e x}}{\sqrt {-c d+b e}}\right )}{\sqrt {c} \sqrt {-c d+b e}}-\frac {(2 b B d-4 A c d+A b e) \text {arctanh}\left (\frac {\sqrt {d+e x}}{\sqrt {d}}\right )}{\sqrt {d}}}{b^3}
\]
[In]
Integrate[((A + B*x)*Sqrt[d + e*x])/(b*x + c*x^2)^2,x]
[Out]
((b*Sqrt[d + e*x]*(b*B*x - A*(b + 2*c*x)))/(x*(b + c*x)) + ((4*A*c^2*d + b^2*B*e - b*c*(2*B*d + 3*A*e))*ArcTan
[(Sqrt[c]*Sqrt[d + e*x])/Sqrt[-(c*d) + b*e]])/(Sqrt[c]*Sqrt[-(c*d) + b*e]) - ((2*b*B*d - 4*A*c*d + A*b*e)*ArcT
anh[Sqrt[d + e*x]/Sqrt[d]])/Sqrt[d])/b^3
Maple [A] (verified)
Time = 0.42 (sec) , antiderivative size = 166, normalized size of antiderivative = 1.05
| | |
method | result | size |
| | |
pseudoelliptic |
\(-\frac {-4 \left (A \,c^{2} d -\frac {3 b \left (A e +\frac {2 B d}{3}\right ) c}{4}+\frac {b^{2} B e}{4}\right ) x \left (c x +b \right ) \sqrt {d}\, \arctan \left (\frac {c \sqrt {e x +d}}{\sqrt {\left (b e -c d \right ) c}}\right )+\left (x \left (c x +b \right ) \left (-4 A c d +b \left (A e +2 B d \right )\right ) \operatorname {arctanh}\left (\frac {\sqrt {e x +d}}{\sqrt {d}}\right )+\sqrt {e x +d}\, \left (\left (-B x +A \right ) b +2 A c x \right ) \sqrt {d}\, b \right ) \sqrt {\left (b e -c d \right ) c}}{\sqrt {\left (b e -c d \right ) c}\, \sqrt {d}\, \left (c x +b \right ) b^{3} x}\) |
\(166\) |
derivativedivides |
\(2 e^{2} \left (-\frac {\frac {\left (\frac {1}{2} A b c e -\frac {1}{2} b^{2} B e \right ) \sqrt {e x +d}}{c \left (e x +d \right )+b e -c d}+\frac {\left (3 A b c e -4 A \,c^{2} d -b^{2} B e +2 B b c d \right ) \arctan \left (\frac {c \sqrt {e x +d}}{\sqrt {\left (b e -c d \right ) c}}\right )}{2 \sqrt {\left (b e -c d \right ) c}}}{e^{2} b^{3}}+\frac {-\frac {A b \sqrt {e x +d}}{2 x}-\frac {\left (A b e -4 A c d +2 B b d \right ) \operatorname {arctanh}\left (\frac {\sqrt {e x +d}}{\sqrt {d}}\right )}{2 \sqrt {d}}}{e^{2} b^{3}}\right )\) |
\(172\) |
default |
\(2 e^{2} \left (-\frac {\frac {\left (\frac {1}{2} A b c e -\frac {1}{2} b^{2} B e \right ) \sqrt {e x +d}}{c \left (e x +d \right )+b e -c d}+\frac {\left (3 A b c e -4 A \,c^{2} d -b^{2} B e +2 B b c d \right ) \arctan \left (\frac {c \sqrt {e x +d}}{\sqrt {\left (b e -c d \right ) c}}\right )}{2 \sqrt {\left (b e -c d \right ) c}}}{e^{2} b^{3}}+\frac {-\frac {A b \sqrt {e x +d}}{2 x}-\frac {\left (A b e -4 A c d +2 B b d \right ) \operatorname {arctanh}\left (\frac {\sqrt {e x +d}}{\sqrt {d}}\right )}{2 \sqrt {d}}}{e^{2} b^{3}}\right )\) |
\(172\) |
risch |
\(-\frac {A \sqrt {e x +d}}{b^{2} x}-\frac {e \left (\frac {\frac {2 \left (\frac {1}{2} A b c e -\frac {1}{2} b^{2} B e \right ) \sqrt {e x +d}}{c \left (e x +d \right )+b e -c d}+\frac {\left (3 A b c e -4 A \,c^{2} d -b^{2} B e +2 B b c d \right ) \arctan \left (\frac {c \sqrt {e x +d}}{\sqrt {\left (b e -c d \right ) c}}\right )}{\sqrt {\left (b e -c d \right ) c}}}{b e}-\frac {\left (-A b e +4 A c d -2 B b d \right ) \operatorname {arctanh}\left (\frac {\sqrt {e x +d}}{\sqrt {d}}\right )}{b e \sqrt {d}}\right )}{b^{2}}\) |
\(175\) |
| | |
|
|
|
[In]
int((B*x+A)*(e*x+d)^(1/2)/(c*x^2+b*x)^2,x,method=_RETURNVERBOSE)
[Out]
-1/((b*e-c*d)*c)^(1/2)/d^(1/2)*(-4*(A*c^2*d-3/4*b*(A*e+2/3*B*d)*c+1/4*b^2*B*e)*x*(c*x+b)*d^(1/2)*arctan(c*(e*x
+d)^(1/2)/((b*e-c*d)*c)^(1/2))+(x*(c*x+b)*(-4*A*c*d+b*(A*e+2*B*d))*arctanh((e*x+d)^(1/2)/d^(1/2))+(e*x+d)^(1/2
)*((-B*x+A)*b+2*A*c*x)*d^(1/2)*b)*((b*e-c*d)*c)^(1/2))/(c*x+b)/b^3/x
Fricas [B] (verification not implemented)
Leaf count of result is larger than twice the leaf count of optimal. 383 vs. \(2 (140) = 280\).
Time = 0.52 (sec) , antiderivative size = 1574, normalized size of antiderivative = 9.96
\[
\int \frac {(A+B x) \sqrt {d+e x}}{\left (b x+c x^2\right )^2} \, dx=\text {Too large to display}
\]
[In]
integrate((B*x+A)*(e*x+d)^(1/2)/(c*x^2+b*x)^2,x, algorithm="fricas")
[Out]
[1/2*(sqrt(c^2*d - b*c*e)*((2*(B*b*c^2 - 2*A*c^3)*d^2 - (B*b^2*c - 3*A*b*c^2)*d*e)*x^2 + (2*(B*b^2*c - 2*A*b*c
^2)*d^2 - (B*b^3 - 3*A*b^2*c)*d*e)*x)*log((c*e*x + 2*c*d - b*e + 2*sqrt(c^2*d - b*c*e)*sqrt(e*x + d))/(c*x + b
)) - ((A*b^2*c^2*e^2 - 2*(B*b*c^3 - 2*A*c^4)*d^2 + (2*B*b^2*c^2 - 5*A*b*c^3)*d*e)*x^2 + (A*b^3*c*e^2 - 2*(B*b^
2*c^2 - 2*A*b*c^3)*d^2 + (2*B*b^3*c - 5*A*b^2*c^2)*d*e)*x)*sqrt(d)*log((e*x - 2*sqrt(e*x + d)*sqrt(d) + 2*d)/x
) - 2*(A*b^2*c^2*d^2 - A*b^3*c*d*e - ((B*b^2*c^2 - 2*A*b*c^3)*d^2 - (B*b^3*c - 2*A*b^2*c^2)*d*e)*x)*sqrt(e*x +
d))/((b^3*c^3*d^2 - b^4*c^2*d*e)*x^2 + (b^4*c^2*d^2 - b^5*c*d*e)*x), -1/2*(2*sqrt(-c^2*d + b*c*e)*((2*(B*b*c^
2 - 2*A*c^3)*d^2 - (B*b^2*c - 3*A*b*c^2)*d*e)*x^2 + (2*(B*b^2*c - 2*A*b*c^2)*d^2 - (B*b^3 - 3*A*b^2*c)*d*e)*x)
*arctan(sqrt(-c^2*d + b*c*e)*sqrt(e*x + d)/(c*e*x + c*d)) + ((A*b^2*c^2*e^2 - 2*(B*b*c^3 - 2*A*c^4)*d^2 + (2*B
*b^2*c^2 - 5*A*b*c^3)*d*e)*x^2 + (A*b^3*c*e^2 - 2*(B*b^2*c^2 - 2*A*b*c^3)*d^2 + (2*B*b^3*c - 5*A*b^2*c^2)*d*e)
*x)*sqrt(d)*log((e*x - 2*sqrt(e*x + d)*sqrt(d) + 2*d)/x) + 2*(A*b^2*c^2*d^2 - A*b^3*c*d*e - ((B*b^2*c^2 - 2*A*
b*c^3)*d^2 - (B*b^3*c - 2*A*b^2*c^2)*d*e)*x)*sqrt(e*x + d))/((b^3*c^3*d^2 - b^4*c^2*d*e)*x^2 + (b^4*c^2*d^2 -
b^5*c*d*e)*x), -1/2*(2*((A*b^2*c^2*e^2 - 2*(B*b*c^3 - 2*A*c^4)*d^2 + (2*B*b^2*c^2 - 5*A*b*c^3)*d*e)*x^2 + (A*b
^3*c*e^2 - 2*(B*b^2*c^2 - 2*A*b*c^3)*d^2 + (2*B*b^3*c - 5*A*b^2*c^2)*d*e)*x)*sqrt(-d)*arctan(sqrt(e*x + d)*sqr
t(-d)/d) - sqrt(c^2*d - b*c*e)*((2*(B*b*c^2 - 2*A*c^3)*d^2 - (B*b^2*c - 3*A*b*c^2)*d*e)*x^2 + (2*(B*b^2*c - 2*
A*b*c^2)*d^2 - (B*b^3 - 3*A*b^2*c)*d*e)*x)*log((c*e*x + 2*c*d - b*e + 2*sqrt(c^2*d - b*c*e)*sqrt(e*x + d))/(c*
x + b)) + 2*(A*b^2*c^2*d^2 - A*b^3*c*d*e - ((B*b^2*c^2 - 2*A*b*c^3)*d^2 - (B*b^3*c - 2*A*b^2*c^2)*d*e)*x)*sqrt
(e*x + d))/((b^3*c^3*d^2 - b^4*c^2*d*e)*x^2 + (b^4*c^2*d^2 - b^5*c*d*e)*x), -(sqrt(-c^2*d + b*c*e)*((2*(B*b*c^
2 - 2*A*c^3)*d^2 - (B*b^2*c - 3*A*b*c^2)*d*e)*x^2 + (2*(B*b^2*c - 2*A*b*c^2)*d^2 - (B*b^3 - 3*A*b^2*c)*d*e)*x)
*arctan(sqrt(-c^2*d + b*c*e)*sqrt(e*x + d)/(c*e*x + c*d)) + ((A*b^2*c^2*e^2 - 2*(B*b*c^3 - 2*A*c^4)*d^2 + (2*B
*b^2*c^2 - 5*A*b*c^3)*d*e)*x^2 + (A*b^3*c*e^2 - 2*(B*b^2*c^2 - 2*A*b*c^3)*d^2 + (2*B*b^3*c - 5*A*b^2*c^2)*d*e)
*x)*sqrt(-d)*arctan(sqrt(e*x + d)*sqrt(-d)/d) + (A*b^2*c^2*d^2 - A*b^3*c*d*e - ((B*b^2*c^2 - 2*A*b*c^3)*d^2 -
(B*b^3*c - 2*A*b^2*c^2)*d*e)*x)*sqrt(e*x + d))/((b^3*c^3*d^2 - b^4*c^2*d*e)*x^2 + (b^4*c^2*d^2 - b^5*c*d*e)*x)
]
Sympy [F(-1)]
Timed out. \[
\int \frac {(A+B x) \sqrt {d+e x}}{\left (b x+c x^2\right )^2} \, dx=\text {Timed out}
\]
[In]
integrate((B*x+A)*(e*x+d)**(1/2)/(c*x**2+b*x)**2,x)
[Out]
Timed out
Maxima [F(-2)]
Exception generated. \[
\int \frac {(A+B x) \sqrt {d+e x}}{\left (b x+c x^2\right )^2} \, dx=\text {Exception raised: ValueError}
\]
[In]
integrate((B*x+A)*(e*x+d)^(1/2)/(c*x^2+b*x)^2,x, algorithm="maxima")
[Out]
Exception raised: ValueError >> Computation failed since Maxima requested additional constraints; using the 'a
ssume' command before evaluation *may* help (example of legal syntax is 'assume(b*e-c*d>0)', see `assume?` for
more detail
Giac [A] (verification not implemented)
none
Time = 0.28 (sec) , antiderivative size = 214, normalized size of antiderivative = 1.35
\[
\int \frac {(A+B x) \sqrt {d+e x}}{\left (b x+c x^2\right )^2} \, dx=-\frac {{\left (2 \, B b c d - 4 \, A c^{2} d - B b^{2} e + 3 \, A b c e\right )} \arctan \left (\frac {\sqrt {e x + d} c}{\sqrt {-c^{2} d + b c e}}\right )}{\sqrt {-c^{2} d + b c e} b^{3}} + \frac {{\left (2 \, B b d - 4 \, A c d + A b e\right )} \arctan \left (\frac {\sqrt {e x + d}}{\sqrt {-d}}\right )}{b^{3} \sqrt {-d}} + \frac {{\left (e x + d\right )}^{\frac {3}{2}} B b e - 2 \, {\left (e x + d\right )}^{\frac {3}{2}} A c e - \sqrt {e x + d} B b d e + 2 \, \sqrt {e x + d} A c d e - \sqrt {e x + d} A b e^{2}}{{\left ({\left (e x + d\right )}^{2} c - 2 \, {\left (e x + d\right )} c d + c d^{2} + {\left (e x + d\right )} b e - b d e\right )} b^{2}}
\]
[In]
integrate((B*x+A)*(e*x+d)^(1/2)/(c*x^2+b*x)^2,x, algorithm="giac")
[Out]
-(2*B*b*c*d - 4*A*c^2*d - B*b^2*e + 3*A*b*c*e)*arctan(sqrt(e*x + d)*c/sqrt(-c^2*d + b*c*e))/(sqrt(-c^2*d + b*c
*e)*b^3) + (2*B*b*d - 4*A*c*d + A*b*e)*arctan(sqrt(e*x + d)/sqrt(-d))/(b^3*sqrt(-d)) + ((e*x + d)^(3/2)*B*b*e
- 2*(e*x + d)^(3/2)*A*c*e - sqrt(e*x + d)*B*b*d*e + 2*sqrt(e*x + d)*A*c*d*e - sqrt(e*x + d)*A*b*e^2)/(((e*x +
d)^2*c - 2*(e*x + d)*c*d + c*d^2 + (e*x + d)*b*e - b*d*e)*b^2)
Mupad [B] (verification not implemented)
Time = 11.03 (sec) , antiderivative size = 2558, normalized size of antiderivative = 16.19
\[
\int \frac {(A+B x) \sqrt {d+e x}}{\left (b x+c x^2\right )^2} \, dx=\text {Too large to display}
\]
[In]
int(((A + B*x)*(d + e*x)^(1/2))/(b*x + c*x^2)^2,x)
[Out]
(atan((((-c*(b*e - c*d))^(1/2)*((2*(d + e*x)^(1/2)*(10*A^2*b^2*c^3*e^4 + 32*A^2*c^5*d^2*e^2 + B^2*b^4*c*e^4 +
8*B^2*b^2*c^3*d^2*e^2 - 6*A*B*b^3*c^2*e^4 - 32*A^2*b*c^4*d*e^3 - 4*B^2*b^3*c^2*d*e^3 - 32*A*B*b*c^4*d^2*e^2 +
24*A*B*b^2*c^3*d*e^3))/b^4 + (((4*A*b^7*c^2*e^4 - 8*A*b^6*c^3*d*e^3 + 4*B*b^7*c^2*d*e^3)/b^6 + ((4*b^7*c^2*e^3
- 8*b^6*c^3*d*e^2)*(-c*(b*e - c*d))^(1/2)*(d + e*x)^(1/2)*(4*A*c^2*d + B*b^2*e - 3*A*b*c*e - 2*B*b*c*d))/(b^4
*(b^3*c^2*d - b^4*c*e)))*(-c*(b*e - c*d))^(1/2)*(4*A*c^2*d + B*b^2*e - 3*A*b*c*e - 2*B*b*c*d))/(2*(b^3*c^2*d -
b^4*c*e)))*(4*A*c^2*d + B*b^2*e - 3*A*b*c*e - 2*B*b*c*d)*1i)/(2*(b^3*c^2*d - b^4*c*e)) + ((-c*(b*e - c*d))^(1
/2)*((2*(d + e*x)^(1/2)*(10*A^2*b^2*c^3*e^4 + 32*A^2*c^5*d^2*e^2 + B^2*b^4*c*e^4 + 8*B^2*b^2*c^3*d^2*e^2 - 6*A
*B*b^3*c^2*e^4 - 32*A^2*b*c^4*d*e^3 - 4*B^2*b^3*c^2*d*e^3 - 32*A*B*b*c^4*d^2*e^2 + 24*A*B*b^2*c^3*d*e^3))/b^4
- (((4*A*b^7*c^2*e^4 - 8*A*b^6*c^3*d*e^3 + 4*B*b^7*c^2*d*e^3)/b^6 - ((4*b^7*c^2*e^3 - 8*b^6*c^3*d*e^2)*(-c*(b*
e - c*d))^(1/2)*(d + e*x)^(1/2)*(4*A*c^2*d + B*b^2*e - 3*A*b*c*e - 2*B*b*c*d))/(b^4*(b^3*c^2*d - b^4*c*e)))*(-
c*(b*e - c*d))^(1/2)*(4*A*c^2*d + B*b^2*e - 3*A*b*c*e - 2*B*b*c*d))/(2*(b^3*c^2*d - b^4*c*e)))*(4*A*c^2*d + B*
b^2*e - 3*A*b*c*e - 2*B*b*c*d)*1i)/(2*(b^3*c^2*d - b^4*c*e)))/((2*(6*A^3*b^2*c^3*e^5 + 32*A^3*c^5*d^2*e^3 - 4*
B^3*b^3*c^2*d^2*e^3 + A*B^2*b^4*c*e^5 - 32*A^3*b*c^4*d*e^4 + 2*B^3*b^4*c*d*e^4 - 5*A^2*B*b^3*c^2*e^5 + 24*A*B^
2*b^2*c^3*d^2*e^3 - 16*A*B^2*b^3*c^2*d*e^4 - 48*A^2*B*b*c^4*d^2*e^3 + 40*A^2*B*b^2*c^3*d*e^4))/b^6 + ((-c*(b*e
- c*d))^(1/2)*((2*(d + e*x)^(1/2)*(10*A^2*b^2*c^3*e^4 + 32*A^2*c^5*d^2*e^2 + B^2*b^4*c*e^4 + 8*B^2*b^2*c^3*d^
2*e^2 - 6*A*B*b^3*c^2*e^4 - 32*A^2*b*c^4*d*e^3 - 4*B^2*b^3*c^2*d*e^3 - 32*A*B*b*c^4*d^2*e^2 + 24*A*B*b^2*c^3*d
*e^3))/b^4 + (((4*A*b^7*c^2*e^4 - 8*A*b^6*c^3*d*e^3 + 4*B*b^7*c^2*d*e^3)/b^6 + ((4*b^7*c^2*e^3 - 8*b^6*c^3*d*e
^2)*(-c*(b*e - c*d))^(1/2)*(d + e*x)^(1/2)*(4*A*c^2*d + B*b^2*e - 3*A*b*c*e - 2*B*b*c*d))/(b^4*(b^3*c^2*d - b^
4*c*e)))*(-c*(b*e - c*d))^(1/2)*(4*A*c^2*d + B*b^2*e - 3*A*b*c*e - 2*B*b*c*d))/(2*(b^3*c^2*d - b^4*c*e)))*(4*A
*c^2*d + B*b^2*e - 3*A*b*c*e - 2*B*b*c*d))/(2*(b^3*c^2*d - b^4*c*e)) - ((-c*(b*e - c*d))^(1/2)*((2*(d + e*x)^(
1/2)*(10*A^2*b^2*c^3*e^4 + 32*A^2*c^5*d^2*e^2 + B^2*b^4*c*e^4 + 8*B^2*b^2*c^3*d^2*e^2 - 6*A*B*b^3*c^2*e^4 - 32
*A^2*b*c^4*d*e^3 - 4*B^2*b^3*c^2*d*e^3 - 32*A*B*b*c^4*d^2*e^2 + 24*A*B*b^2*c^3*d*e^3))/b^4 - (((4*A*b^7*c^2*e^
4 - 8*A*b^6*c^3*d*e^3 + 4*B*b^7*c^2*d*e^3)/b^6 - ((4*b^7*c^2*e^3 - 8*b^6*c^3*d*e^2)*(-c*(b*e - c*d))^(1/2)*(d
+ e*x)^(1/2)*(4*A*c^2*d + B*b^2*e - 3*A*b*c*e - 2*B*b*c*d))/(b^4*(b^3*c^2*d - b^4*c*e)))*(-c*(b*e - c*d))^(1/2
)*(4*A*c^2*d + B*b^2*e - 3*A*b*c*e - 2*B*b*c*d))/(2*(b^3*c^2*d - b^4*c*e)))*(4*A*c^2*d + B*b^2*e - 3*A*b*c*e -
2*B*b*c*d))/(2*(b^3*c^2*d - b^4*c*e))))*(-c*(b*e - c*d))^(1/2)*(4*A*c^2*d + B*b^2*e - 3*A*b*c*e - 2*B*b*c*d)*
1i)/(b^3*c^2*d - b^4*c*e) - (((d + e*x)^(1/2)*(A*b*e^2 - 2*A*c*d*e + B*b*d*e))/b^2 + ((2*A*c*e - B*b*e)*(d + e
*x)^(3/2))/b^2)/((b*e - 2*c*d)*(d + e*x) + c*(d + e*x)^2 + c*d^2 - b*d*e) - (atanh((4*B^3*c*d^(1/2)*e^4*(d + e
*x)^(1/2))/(2*A*B^2*c*e^5 - (8*A^3*c^3*e^5)/b^2 + 4*B^3*c*d*e^4 + (2*A^3*c^2*e^6)/(b*d) - (16*A*B^2*c^2*d*e^4)
/b + (16*A^2*B*c^3*d*e^4)/b^2) + (2*A^3*c^2*e^6*(d + e*x)^(1/2))/(d^(3/2)*((2*A^3*c^2*e^6)/d - (8*A^3*c^3*e^5)
/b - 16*A*B^2*c^2*d*e^4 + 2*A*B^2*b*c*e^5 + 4*B^3*b*c*d*e^4 + (16*A^2*B*c^3*d*e^4)/b)) - (8*A^3*c^3*e^5*(d + e
*x)^(1/2))/(d^(1/2)*(2*A*B^2*b^2*c*e^5 - 8*A^3*c^3*e^5 + 16*A^2*B*c^3*d*e^4 + 4*B^3*b^2*c*d*e^4 + (2*A^3*b*c^2
*e^6)/d - 16*A*B^2*b*c^2*d*e^4)) - (16*A*B^2*c^2*d^(1/2)*e^4*(d + e*x)^(1/2))/((2*A^3*c^2*e^6)/d - (8*A^3*c^3*
e^5)/b - 16*A*B^2*c^2*d*e^4 + 2*A*B^2*b*c*e^5 + 4*B^3*b*c*d*e^4 + (16*A^2*B*c^3*d*e^4)/b) + (16*A^2*B*c^3*d^(1
/2)*e^4*(d + e*x)^(1/2))/(2*A*B^2*b^2*c*e^5 - 8*A^3*c^3*e^5 + 16*A^2*B*c^3*d*e^4 + 4*B^3*b^2*c*d*e^4 + (2*A^3*
b*c^2*e^6)/d - 16*A*B^2*b*c^2*d*e^4) + (2*A*B^2*c*e^5*(d + e*x)^(1/2))/(d^(1/2)*(2*A*B^2*c*e^5 - (8*A^3*c^3*e^
5)/b^2 + 4*B^3*c*d*e^4 + (2*A^3*c^2*e^6)/(b*d) - (16*A*B^2*c^2*d*e^4)/b + (16*A^2*B*c^3*d*e^4)/b^2)))*(A*b*e -
4*A*c*d + 2*B*b*d))/(b^3*d^(1/2))